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1 Introduction 



The computation of higher order corrections to multi-scale processes is plagued by the presence 
of unstable particles in loop integrals. Formally, a clean description would require a Dyson resum- 
mation of self-energy insertions in order to preserve unitarity [1]; in the context of the Standard 
Model and its extensions, however, the consequent mixing of perturbative orders clearly compro- 
mises gauge invariance. In addition to the unitarity issue, a practical problem is represented by the 
appearance of unphysical threshold singularities at the amplitude level for physical observables. 

In this letter we focus on two-loop electroweak corrections to Standard Model Higgs production 
through gluon fusion, gg H, and decay into two photons, — > 77 (details of the calculation will 
be given in a forthcoming paper [2] ) . We discuss the singular behavior of the amplitudes around the 
normal thresholds induced by internal unstable particles, directly related to two-particle unitarity 
cuts. Concerning the ^ 77 case, note that threshold effects have been analyzed in Ref. [3], with 
special emphasis on the presence of bound states, which are not the subject of our study. 

The singular behavior can be cured trading the real masses for unstable particles, used as 
experimental input data, with the associated complex poles. As shown by the authors of Ref. [4], 
the replacement has to be performed also at the level of the couplings, leading to the so-called 
complex-mass scheme. Consequently, the one- and two-loop integrals needed for the computation 
have to be evaluated with complex internal arguments. 

A minimal implementation of the complex-mass scheme for two-loop electroweak corrections to 
the — > 77 process has been realized in Ref. [5]. The two-loop amplitude is splitted in gauge- 
invariant divergent and non-divergent terms, and complex masses are introduced in the divergent 
part. This solution is formally satisfactory, since the amplitude is finite also at thresholds. Further- 
more, as a consequence of the cancellation mechanism taking place among divergent contributions, 
this scheme does not require any analytic continuation of two-point functions connected to mass 
renormalization to the second Riemann sheet. 

As a drawback, artificially large numerical effects arise around normal thresholds. The issue is 
relevant when we consider Higgs production through gluon fusion at the LHC, where a measurement 
of the Higgs mass can be performed at the per-mille level [6]. The statement that next-to- leading 
order electroweak corrections are known [7] does not appear to be fully adequate, given the possible 
occurrence of large threshold effects of about 10 % with respect to the leading order result. 

In this letter we show that threshold effects for gg ^ H and — > 77 are well under control 
in our improved calculational scheme [2], where a complete implementation of the complex-mass 
scheme at two loops is performed along the lines of Ref. [8]. For the analysis in hadron-hadron 
collisions we refer to Ref. [9]. 

2 Radiative corrections with unstable particles 

In this section we briefly summarize aspects of selected solutions for dealing with unstable 
particles in tree and one-loop calculations; next, we analyze the salient features of a two-loop 
computation. In Section [3] we will study where the numerical impact of complex masses is most 
relevant at the two-loop level by looking into the presence of unphysical infinities and cusps in 
two-loop amplitudes. 

At tree level, if the external legs of a given amplitude are divided into two disjoint subsets and 
if the total quantum numbers of each subset allow for the exchange of a known particle of mass M, 
then the amplitude has a pole satisfying = M^, where p is the four- momentum of the exchanged 
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particle. In leading order (LO) calculations, the masses of these particles, like W and Z bosons, 
are replaced by the location of the poles in the complex p'^ plane. However, the principle of gauge 
invariance must not be violated: in particular. Ward- Slavnov- Taylor (hereafter WST) identities [10] 
have to be preserved, otherwise theoretical uncertainties may get out of control. The incorporation 
of finite-width effects in the theoretical predictions for Lep2 processes was a typical example and, 
at the time, it was argued that the preferable scheme consists in resummation of fermion one-loop 
corrections to vector-boson propagators and inclusion of all remaining fermion one-loop corrections, 
in particular those to the Yang- Mills vertices [11]. 

A possible solution at the next-to-leading order (NLO) level consists in replacing everywhere 
the squared real masses (m^) with the complex poles (sp), couplings included; this is known in 
the literature as complex-mass scheme [4]. Since WST identities arc algebraic relations satisfied 
separately by real and imaginary parts, one starts from WST identities with real masses, satisfied 
at any given order, and replaces everywhere m? — > Sp without violating the invariance. 

In tTirns, this scheme violates unitarity: one cannot identify the two sides of any cut diagram 
with T and respectively (the transition matrix T is defined in terms of the S matrix as S" = 
1 + iT). To summarize, the analytical structure of the S matrix is correctly reproduced when we 
use propagator factors s — Sp, where s is a generic invariant, but unitarity of S requires more, a 
dressed propagator [8,12]. However, we expect that unitarity- violating terms are of higher order; 
in principle, the violating terms should not be enhanced because WST identities are preserved. 

Another drawback of the scheme is that all propagators for unstable particles will have the same 
functional form both in the time-like and space-like regions, while, for a dressed propagator, the 
presence of a pole on the second Riemann sheet does not change the real character of the function 
if we are in a t channel. 

Typical examples of one-loop calculations that require the introduction of complex poles are 
those for processes, like e'^e" — 4f, where part of the amplitude (the so-called signal) factorizes 
into production (g) decay of one or more particles; or processes involving off-shell W -pair production, 
leading to the so-called Coulomb singularity. Another approach in this context is represented by 
the use of effective field theory methods (see recent applications in Ref. [13]). 

In some sense the complex-mass scheme becomes more appealing when we go beyond one loop, 
as described in detail in Ref. [8], where a recipe was designed to derive loop amplitudes out of a 
skeleton expansion. Let m? be the squared bare mass, Sp the corresponding complex pole and S 
the corresponding self-energy: to the requested order we replace everywhere m? with Sp + S(sp), 
which is real by construction. If only one loop is needed, then Sp everywhere (therefore 

justifying the name complex mass). 

Note that the on-shell mass is related to the zero of the real part of the inverse propagator; 
beyond one loop this would show a clash with gauge invariance, since only the complex poles do 
not depend on gauge parameters to all orders. As a consequence, renormalization equations change 
their structure. 

Furthermore, there is also a change of perspective with respect to one-loop calculations. There 
one considered the on-shell masses as input parameters independent of complex poles and derived 
the latter in terms of the former [11]. Here the situation changes: renormalization equations are 
written for real renormalized parameters and solved in terms of (among other things) experimental 
complex poles; the latter have to be reconstructed from on-shell pscudo-observablcs. 

Having described the general setup for a gauge-invariant formulation of unstable particles at 
the multi-loop level, we want to understand where, in a two-loop calculation, the numerical impact 
of complex masses is most relevant. 
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3 Two-loop amplitudes and normal thresholds 

In this section we explore the singular behavior of massive two-loop amplitudes around normal 
thresholds, investigating the origin of square-root and logarithmic singularities. It is worth noting 
that pseudo-thresholds are always outside the physical region. 

Normal thresholds are directly related to unitarity cuts as illustrated in Fig. [H and correspond 
to the leading Landau singularities [14] of self-energy diagrams. When diagrams with more than 
two external legs are present, normal thresholds show up as sub-leading singularities; this can be 
easily understood observing that all diagrams in Fig. [1] are reduced to self-energy configurations 
after shrinking a line which does not intersect any cut to a point. 




Figure 1: Sample cut diagrams for — > 77 showing normal thresholds for = 2M^, 2M^, 2Mf. 

An interesting question is the role played by the leading singularity (the so-called anomalous 
threshold) [15]: already at one loop, there are cases where a non-integrable singularity associated 
with the leading Landau singularity requires the introduction of complex masses [16]. Concerning 
the processes — > 77 and gg — > H, we have verified that the leading Landau singularities of all 
two-loop vertex diagrams lie outside the physical region, and we will drop this issue in the following. 

For our discussion, it is useful to decompose the — > 77 amplitude as 

^=(V2C,Mj)"^^(.4- + m.4-o), (1) 

where Gp is the Fermi-coupling constant, is the mass of the W boson and a is the fine- 
structure constant; and ^^^^ denote the leading order and next-to-leading order amplitudes. 
Furthermore, ^^'^^ can be written as 

where A"^^ is given by the sum of all pure two-loop diagrams, _4^^^ follows from the inclusion of 
the one-loop Higgs wave-function renormalization (WFR) factor and ^^^'^ stems from one-loop 
renormalization of the masses and the Fermi coupling. Note that the U(l) Ward identity forces 
the electromagnetic coupling to go unrenormalized once the external on-shell photons are provided 
with their WFR factors. For electroweak corrections, the amplitude for gg ^ H is given by Eq.([T]), 
with a replaced by the strong-coupling constant ag(/i^) at the renormalization scale /U^. 

3.1 Square-root singularities 

A square-root singularity is represented by a term containing a single inverse power of the 
threshold factor 

(3i = yJl-AMf/Ml, with Mi = M^,M^,Mt. (3) 
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In Ref. [5] it has been shown that square-root singularities are related to: i) derivatives of 
two-point one-loop functions, associated with Higgs WFR; ii) derivatives of three-point one-loop 
functions, generated by mass renormalization; iii) genuine irreducible two-loop diagrams containing 
a one-loop self-energy insertion. Therefore, in general, all three terms of Eq.dJ]) can show a 
behavior. 

Concerning the Higgs WFR factor at one loop, we deal with the four mass patterns shown in 
Fig. [21 note that tadpole diagrams do not affect the threshold behavior. For the top-quark diagram, 
the coefficient of the derivative of the two-point one-loop function contains a positive power of the 
threshold factor Pt; in other words, this diagram is /3f-protected at threshold. Consequently, ^wfr 
of Eq.dJ]) contains, for both processes ^ 77 and gg H, square-root singularities only at the 
2 and 2 thresholds. 




W Z H 



Figure 2: Sample diagrams for the Higgs wave-function renormalization factor at one loop. 

We consider now genuine two-loop diagrams containing a self-energy insertion; they natu- 
rally join terms induced by one-loop mass renormalization as shown in Fig. [3l where bosonic and 
fermionic diagrams are illustrated. The bosonic component is obviously peculiar of the if — > 77 
decay; in addition, we observe that only charged bosonic diagrams are present, because of the 
nature of the triple non-abelian gauge coupling in the Standard Model. 




Figure 3: Two- loop and mass-renormalization diagrams relevant for the analysis of square-root 
singularities. Gray circles represent the sum of all one-loop two-point diagrams; black dots denote 
a derivative. 

Fermionic diagrams are /?t-protected at threshold, and do not require any special care. The 
two-loop vertex containing a W self-energy insertion, and the associated PF-mass renormalization 
contribution, lead instead to a /3,^^-divergent behavior. However, the two-loop irreducible diagram 
of Fig. [3^ can be cast in a representation where the singular part is completely written in terms of 
the one-loop diagrams of Fig. [3b- Moreover, it is possible to check explicitly that the unphysical 
behavior, generated by the two-loop diagram of Fig. [3^, exactly cancels the divergency 
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due to one-loop Vl^-mass renormalization of Fig. [8)3 (performed in the complex-mass setup we are 
going to describe in Section H]). 

This cancellation mechanism corroborates the general picture of Ref. [1]: self-energy insertions 
signal the presence of an unstable particle, and are the consequence of a misleading organization of 
the perturbative expansion; Dyson-resummed propagators should be used and complex poles should 
replace real masses as input data. The outcome of our analysis, performed at the algebraic level, 
is straightforward: at the amplitude level, only the inclusion of the Higgs WFR factor generates 
square-root singularities at the 2M^^ and 2M^ thresholds for both processes — > 77 and gg — > H. 

3.2 Logarithmic singularities 

We briefly address here the issue of logarithmic singularities, given by terms containing a factor 
ln{—Pf — zO), generated by the diagrams of Fig. [H As thoroughly discussed in Ref. [5], the scalar 
configuration associated with the diagrams illustrated in Fig. U] generates a logarithmic singularity. 
If the massive loop is made of top quarks, the scalar integral appears at the amplitude level with 
a multiplicative factor and the logarithmic singularity is /3|-protected at threshold. The same 
consideration is not valid for a W loop; here the logarithmic singularity can be viewed as the 
remnant of a Coulomb singularity in the one-loop sub-diagram. 




'o'o'a'a'O"- g 




Figure 4: Irreducible two-loop vertex diagrams which can generate a logarithmic divergency. 

4 Complex masses 

A pragmatic gauge-invariant solution to the problem of threshold singularities due to unstable 
particles for the 77 ^ 77 decay has been introduced and formalized in Ref. [5] (ad hoc introduction 
of a width is however common practice in the literature). After reviewing the corresponding scheme, 
termed in the following as minimal complex-mass (MCM) scheme, we will discuss the extension 
to the full complex-mass (CM) scheme, looking ahead to precise predictions for the production 
mechanism gg ^ H. 

4.1 Minimal complex-mass scheme 

In the MCM scheme the NLO amplitude of Eq.([l]) and Eq.([2]) is decomposed according to 

^NLO ^ ^ ^ + A^^^ In (_^2^ _ ,0) + AnEU, (4) 
i=W,Z 

where square-root- (AsR,j ) and logarithmic-singular (j4log) terms have been isolated from the 
component which is finite for /3j — > (^rem)- From the discussion of Section [3l we know that 
^SR,i, with i = W,Z, is generated for both — > 77 and gg ^ H hy Higgs WFR at one loop; j4log 
shows up for — > 77 only, and is induced by the bosonic diagram of Fig. [H 
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After proving that all coefficients in Eq.Q, gauge-parameter independent by construction, 
satisfy the WST identities, we minimally modify the amplitude introducing the complex-mass 
scheme of Ref. [4] for the divergent terms. In principle, two steps are required: first, the real 
masses of the W and Z bosons, used as input data, are replaced by the corresponding complex 
poles in the threshold factors Pi, i = W,Z, and in the coefficients A^R^i and ^log (also at the 
level of the couplings); second, the real parts of the W and Z self-energies stemming from mass 
renormalization at one loop are traded for the complete self-energies, including imaginary parts. 

In practice, the second step amounts to a replacement of the conventional on-shell mass renor- 
malization equation with the associated expression for the complex poles, 



Mf 



ReSf)(M2) 



(1) 



(5) 



where t!i\m1), with i = W,Z, denotes the W (Z) one-loop self-energy, and complex poles are 
defined as 



^ii il^i - in) , ^4 = ^i 



1 



2Mt 



(6) 



Here Mj and Fj are the canonical on-shell values for the mass and the width of unstable gauge 
bosons. 

Note that, concerning the W boson, the replacement of the real part of the self-energy with the 
full expression has to be performed also at the level of the Fermi-coupling renormalization equation, 
which becomes 



1/2 



^ GpSyy ^ 



4^/2 



2si 



Incl (7) 



Here g is the bare (or MS'-renormalized) weak-coupling constant, and the squared cosine of the 



weak-mixing angle Cg = cos d {sg = sin 6 = 1 — Cg) is fixed by Cg 



Here it is important to note that the cancellation mechanism between two-loop diagrams and 
one-loop mass-renormalization terms generated by the decomposition of Fig. [3^,b, mentioned com- 
menting Eq.Q, is a consequence of the introduction of Eq.Q and does not take place in the 
conventional on-shell renormalization framework for real masses. This cancellation has a striking 
consequence: one-loop mass renormalization contributes only to the remainder ^rem; and the steps 
summarized in Eq.([5]) and Eq.([7I) are not required anymore. 

The MCM scheme allows for a straightforward removal of unphysical infinities: real masses of 
unstable gauge bosons are traded for complex poles in divergent terms, gauge-parameter invariance 
and WST identities are preserved and the amplitude has a decent threshold behavior, as shown in 
Fig. [5] for the NLO electroweak corrections to the i7 ^ 77 decay width. Here the dotted curve 
represents the result obtained using conventional on-shell masses for unstable gauge bosons as input 
data; the two-loop amplitude is artificially infinite at threshold and badly fails to approximate the 
correct result above threshold, as a consequence of the severe behavior of Eq.(|4|) (enhanced 
above threshold by the fact that the compensation illustrated in Fig. [3^,b does not occur in the 
pure real-mass setup). 

A nice feature of the MCM scheme (dashed curve) is its simplicity: as a consequence of the 
cancellation mechanism taking place among divergent contributions, this scheme does not require 
any analytic continuation of two-point functions connected to mass renormalization to the unphys- 
ical Riemann sheet, because the replacement indicated in Eq.Q is not needed. The MCM scheme, 
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Figure 5: Percentage NLO electroweak corrections to the partial width for H — > 77; here r ^ = 
r^°(l + 5Ew)- MCM (CM) scheme is described in SectionO(4.2). Setup is described in Section[5l 



however, does not deal with cusps associated with the crossing of normal thresholds, as shown in 
Fig. [5] for the WW threshold. 

4.2 Complex-mass scheme 

The large and artificial effects arising around normal thresholds in the MCM scheme or in 
a scheme where the masses of unstable particles are kept real (in this case the amplitude exactly 
diverges at threshold) are aesthetically unattractive. In addition, they represent a concrete problem 
in assessing the impact of two-loop electroweak corrections on processes relevant for the LHC. An 
important example is represented by Higgs production via gluon fusion: here, the large effect of 
NNLO QCD corrections naturally suggests to investigate possible sources of additional sizeable 
corrections. 

We have therefore undertaken the task of introducing the complete complex-mass scheme of 
Ref. [4], as explained for a two-loop calculation in Ref. [8] (see also Ref. [12]), for evaluating 
two-loop electroweak corrections to the Standard Model processes if — > 77 and gg ^ H. The 
procedure described in Section 14.11 for the divergent terms of Eq.(j4]) has been extended to the 
remainder Arem- In particular, all two- loop diagrams have been computed with complex masses 
for the internal vector bosons relying on the techniques developed in Ref. [2, 17]. 

In the full CM setup, the real parts of the W and Z self-energies induced by one-loop renormal- 
ization of the masses and the couplings have to be traded for the associated complex expressions 
by means of Eq.Q and Eq.Q. However, we notice that: 

i) for H — > 77, the Z boson self-energy connected to mass renormalization would arise only 
from the tree-level couplings of the photons, entailing an overall factor g'^Sg for the LO amplitude; 
in the CM scheme, in fact, Sg is expressed through the ratio of the vector-boson masses. However, 
because of our choice for the input-parameter set [5], the factor g'^Sg is re-absorbed by introducing 
the fine-structure constant a. In addition, as we said earlier, the electromagnetic coupling goes 
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unrenormalized, once WFR factors for on-shell photons are included. 

ii) for gg ^ H, Z mass renormalization clearly does not play any role. In addition, the tree- 
level coupling of the Higgs field to the top quark contains a factor g/M^; after combining the 
renormalization of the weak-coupling constant g, related to the Fermi-coupling constant through 
Eq.©, with mass renormalization for the W boson, encoded in Eq.([5]), the W self-energy evaluated 
at the complex pole drops out (see also Ref. [18]). 

As a result, for gg ^ H it is enough to replace the real masses for the W and Z bosons with 
their complex poles, as usual also in the couplings; for H 77, one has also to trade the real part of 
the W self-energy for its full complex expression at the level of mass and coupling renormalization, 
via Eq.(l5|) and Eq.©. The effect for the ^ 77 decay mode, shown in Fig. [5] (solid curve), is 
a full smoothing of the unphysical cusp associated with the WW threshold; although numerical 
negligible, it provides a benchmark for the gg ^ H study we will perform in Section [5l 

The scheme can be easily extended to the fermionic sector replacing also the top-quark real 
mass by its complex pole. Since the behavior associated with the ti threshold is f3t -protected, we 
do not purse this issue here. 

The introduction of gauge-invariant complex poles for gauge bosons leads to technical com- 
plications, due to the fact that one- and two-loop integrals have to be computed with complex 
arguments. If only internal masses are complexified, the analytical continuation of loop integrals 
does not pose any additional problem; after writing the parametric representation of one- and two- 
loop diagrams, it is easy to control that squared masses have semi-positive definite coefficients; 
therefore the replacement — iO^Sp = fj? — ifi'j is straightforward. 

One-loop two-point functions arising in the reduction of the amplitude, instead, have to be 
carefully treated; here the external squared momentum can be complex and logarithms have to be 
extended to the second Riemann sheet. In general, the presence of complex momenta in two- loop 
diagrams demands an analytical continuation also for poly logarithms. In all cases, the correct 
analytical continuation is determined by the request that the value for a stable gauge boson should 
be smoothly approached when the coupling tends to zero. This is achieved starting from a complex 
argument, z = Zji + izj, defining z = Zji — iO, and replacing ordinary logarithms and polylogarithms 
with 

In (z ; z) = In z — 2 i vr ("-^s) , Li2 (z ; z) = Li2 (z) — 2 in In z^ 9 {z^ — 1) . (8) 

4.3 External unstable particles 

From a formal perspective, external unstable particles should not appear in any computation, 
since they cannot be included in the asymptotic states forming the bases of the Hilbert space. 
Concerning Higgs physics, however, available calculations deal with an external on-shell Higgs 
boson and do not perform the ultimate step of introducing a complex pole and the associated 
residue for the decaying Higgs [19]. 

For the production process gg H, we have verified that there are no practical problems 
associated with gauge-parameter invariance and WST identities once we deal with an on-shell Higgs. 
For the — > 77 decay, there is a LO contribution containing the bare Higgs mass, represented 
by a charged Higgs-Kibble one-loop triangle diagram. Standard mass renormalization introduces 
the on-shell Higgs mass, through the real part of the one-loop two-point self-energy, and leads to a 
violation of the WST identities above the WW threshold. In both aforementioned MCM and CM 
schemes the real part of the Higgs self-energy stemming from mass renormalization is traded for 
the complex expression, even if the external Higgs boson is assumed to be an on-shell particle. 
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5 Numerical effects for Higgs production 



The dominant production mechanism of the Standard Model Higgs boson at the LHC is the 
gluon fusion process, gg — > H: in this section we discuss the numerical impact at the two-particle 
vector-boson thresholds of the two-loop electroweak corrections in the framework of the MCM and 
CM schemes. 

In the computation we have set light-fermion masses to zero and introduced the W and Z boson 
complex poles by means of Eq.®. As input parameters we have used the following values taken 
from Ref. [20]: 

= 80.398 GeV, = 91.1876 GeV, 

Tz = 2.4952 GeV, Gp = 1.16637 x 10"^ QeV-^. 

For the mass of the top quark, we have used Mt = 170.9 GeV [21]; for the width of the W boson, 
we have chosen the value F^y = 2.093 GeV, predicted by the Standard Model with electroweak and 
QCD corrections at one loop. 

Our results for (^ewj defined through fj^^ = a^^{l + 6Ew), are shown in Fig. [U where we include 
the complete corrections, comprehensive of light- and top-quark contributions, comparing MCM 
and CM scheme. The corresponding numerical results are given in Tab. [TJ 




In Fig. [7] we show the details of the region around the WW threshold, including the result 
obtained using purely real masses. 

The numerical relevance of NLO electroweak corrections around the vector-boson thresholds 
depends crucially on the implementation of the renormalization scheme: the relative corrections 
in the MCM scheme reach about 10% (2%) at the WW {ZZ) threshold; in the CM scheme, they 
amount to 5% at the WW threshold and vanish at the ZZ one. 
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Table 1: Percentage electroweak corrections in MCM (Section 14. 1|) and CM (Section 14. 2|) schemes 
as a function of the Higgs mass in GeV {Alt = 170.9 GeV). 
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Figure 7: Same as Fig.[6l around the WW threshold. MCM (CM) scheme is described in Section HTTl 
(4.2). Setup of Section [5l 



6 Conclusions 



In this paper we have considered the extension of the complex-mass scheme to two-loop multi- 
scale calculations. After discussing the general setup we have given numerical results for the 
two-loop percentage corrections to the i7 — > 77 and gg ^ H processes in the Standard Model 
around the vector-boson thresholds. 

We have compared a minimal implementation of the complex-mass scheme and the complete 
one. For Higgs masses close to the WW and ZZ thresholds, NLO electroweak corrections can be 
considered under control only after a full implementation of the complex-mass scheme. 

The electroweak scaling factor for the cross section does not exceed the +6% level in the 
range 100 GeV < M^ < 200 GeV; incongruent +10% effects around thresholds are avoided, as a 
consequence of the complex-mass scheme employed. 
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